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Abstract 

The Global Colour Model (GCM) of QCD is a quark-gluon quantum field 
theory that very successfully models QCD for low energy hadronic processes. 
An effective gluon correlator models the interaction between quark currents. 
Functional integral calculus allows the GCM to be hadronised (i.e. expressed in 
terms of meson and baryon fields). The dominant configuration of the hadronic 
functional integrals is revealed to be the constituent quark effect, and is identi- 
cal to one version of the truncated quark Dyson-Schwinger equations (tDSE). 
However the GCM shows that hadronic physics requires processes that go be- 
yond the tDSE. In this review examples of meson and nucleon processes are 
given. The GCM also plays a pivotal role in showing how QCD may be related 
to many other hadronic models 
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1 Introduction 

We review the Global Colour Model (GCM) (Cahill and Roberts (1985) §) of Quantum Chro- 
modynamics (QCD) with particular emphasis on its hadronisation and the resulting applications 
to low energy meson and nucleon processes. Other reviews of the GCM are Cahill (1989) H, 
Cahill (1992)§, Roberts and Williams (1994)j§, Tandy (1996)@, Cahill and Gunner (1997)@ 
and Tandy (1997)0. QCD is defined by the quantisation of the quark-gluon fields with 'classical' 
action S\q, q, A a ^\. However all evidence for the quarks and gluons is provided by the properties and 
interactions of the hadrons and by processes involving the electroweak particles. These hadronic 
laws are encoded in an effective action S[ir, p, cj, .., N, N, .], where ir(x), .., N(x), N(x), ... are fields 
describing composite constituent (equivalcntly core or bare) hadrons, with 'centre-of-mass coordi- 
nates' x. These hadronic fields are to be quantised subject to this effective action, yielding finally 
the observables of QCD. Such an effective action must clearly be non-local because of the com- 
posite nature of these hadrons. While a general derivation of S[tt, .., N, N, .} from S[q, q, A^] has 
not been achieved, it is possible to do this hadronisation within the GCM. The hadronisation uses 
Functional Integral Calculus (FIC) techniques, which amount to dynamically determined changes 
of integration variables in the functional integral formulation of the GCM; 

J DqDqDAexp(SGCM[A,q,q])= J Dir..DNDN.. exp(-S had [TT, .., N, N, ..}). (1) 

Here the constituent hadrons are essentially the normal modes. A particular feature of the GCM 
is that it plays a pivotal role in relating various seemingly different modellings of QCD as shown in 
Fig.0, and these relationships will be discussed later. A key task in using the GCM is the determi- 
nation of the low energy quark-gluon processes from experimental data. In recent years there have 
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been three such extractions of increasing complexity: GCM95§, GCM97§ and GCM98@. We 
report here new detailed ab initio studies of the nucleon within the GCM in which one proceeds 
in a rigorous manner from the experimentally determined low energy quark-gluon processes to 
detailed dynamical studies of the nucleon, including dressing by the pions. 
Truncated SDE »- 
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Figure 1: Relational map of the GCM to QCD and various other modellings including the Nambu - Jona- 
Lasinio (NJL), Quantum Hadrodynamics (QHD), Quark Meson Coupling (QMC) and Chiral Perturbation 
Theory (ChPT). 



2 Global Colour Model 

Here we discuss the construction of the GCM from QCD. In the functional integral approach 
correlators are defined by 

_ J DqDqDADCDC^.qjx) exp(-S QC£> [Ag,g, 

J DqDqDADCDCcxp(-S Q cD[A,q : q,C,C}) [ ' 

where the 'classical' action defining chromodynamics is, in Euclidean metric, 
S QCD [q,q,A«] = J d 4 x(^F* u F« u + ^{d^f 

+ q(^(d»-ig^A«)+M) q y (3) 

This involves the quark and gluon fields and the field strength tensor for the gluon fields. A4 = 
{m u , rrid, ■ ■} are the quark current masses, and ghost (C, C) and gauge fixing terms must be added 
to Sqcd in (§)■ The chromodynamic action clearly has two important invariance groups, Poincare 
symmetry and the local colour symmetry. The various complete (denoted by scripted symbols) 
correlators Q lead to experimental observables. They are related by an infinite set of coupled 
Dyson-Schwinger Equations (DSE), and by the Slavnov- Taylor gauge-symmetry identities and, 
in the chiral limit, to the axial Ward-Takahashi identity (AWTI). The usual truncation of these 
DSE causes the violation of all these identities. The correlators in (||) may be extracted from the 
generating functional of QCD, 

Zqcd[v, V, J\ = I DqDqDADCDCexp(-S QC D[A,q, q,C, C]+fjq + QV+ J A). (4) 

Functional transformations and approximations lead to the GCM; briefly and not showing source 
terms for convenience, the gluon and ghost integrations are formally performed 

J DqDqD ADC X>Ccxp(- S QC d) = J DqDqcxp(- J q(j.d + M)q 

+f / 3 a ,(xW(y)G,Ax ~y) + §J fflifiGfi, + ), (5) 
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where j^(x) — q{x)^^ fi q{x), go is the bare coupling constant, and Q liU {x) is the gluon correlator 
with no quark loops but including ghosts (C, C) 



G^(x - y) 



J DAVCDCA a Jx)Al{y)eMSQCD[A, C, C]) 



(6) 



/ DADCDCexp(-S QC D[A, C, C]) 

Fig.|| shows successive terms in d). This infinite sequence is a direct consequence of the local 
non-abelian colour symmetry. The terms of higher order than the term quartic in the quark 
fields arc difficult to explicitly retain in any analysis. The GCM models the effect of higher order 
terms by replacing the coupling constant go by a momentum dependent quark-gluon coupling g(s), 
and neglecting terms like Q^ p and higher order in (||). This g(s) is a restricted form of vertex 
function. The modification goGfiv(p) —> D^fjj) — g{p 2 )Q \iv(jp)g{p 2 ) and the truncation then defines 
the GCM. We also call this effective quark-quark coupling correlator D llv {p) the effective gluon 
correlator and show included processes in Fig|^. As discussed here D iLV {p) may be determined 
from experimental data. We then formally define the GCM as the quark-gluon field theory with 
the action 



Sgcm[<1, q, A*] 



qfr.d + M 



iA a ^ Jfl )q+UlD-Z(id)A a v 



and the generating functional 



Z[J,f),v} = J DqDqDAexp(-SaCM[q,q,A^]+f)q + qr] + J^Al). 



(7) 



(8) 






(a) (b) (c) 

Figure 2: Diagrammatic representation of successive terms in the quark action in (ISf) - The quark-gluon 
vertex strength is go, while the gluon-gluon vertices (including gluon correlators) are fully dressed except 
for quark loops. 





(a) 



(b) 



(c) 



Figure 3: (a) The GCM effective _D M „ in (J7j) , (b) example of correlations formally included in (a), and in 
(c) an n = 4 process not formally included in (a), but which is modelled in the GCM via the specific form 
of D av ■ 



Here D^{p) is the matrix inverse of the Fourier transform of D^ u (x). The action Sgcm is 
invariant under q — > Uq, q — > qW and A^X a — > UA®\ a W, where U is a global 3x3 unitary 
colour matrix; this is the global colour symmetry of the GCM. The gluon self-interactions that 
arise as a consequence of the local colour symmetry in (||) and the ghost and vertex effects lead 
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to D~^(p) in (0) being non-quadratic. Hence, in effect, the GCM models the QCD local gluonic 
action J F" [A] F^ v [A] in Sqcd of @) which has local colour symmetry, by a highly nonlocal action 
in the last term of (0) which has global colour symmetry. It is important to appreciate that while 
the GCM has a formal global colour symmetry, the detailed dynamical consequences of the local 
colour symmetry of QCD are modelled by the particular form of D(s). There is an Infrared (IR) 
saturation effect which, in conjunction with the dynamical breaking of chiral symmetry, appears 
to suppress details of the formal colour gauge symmetry of QCD. As well, in the chiral limit 
A4 — > 0, the GCM action has Ul(Nf) (g> Ur(Nf) symmetry: the q — q part of S may be written 

= QrI^QR + Ql7i^Ql, where qn : L = Pnxq and q RL = qPb,R- These two parts are separately 
invariant under q R — > Unqn, q R — > q R U R and q L — > C/l<?l, ~q~L ~ * Il^l- 

A key feature of the GCM analysis is the demonstration that this effective gluon correlator 
Duvijp) is successful for a variety of hadronic processes, i.e. it is a universal feature of low energy 
hadronic processes. The success of the GCM appears to be based on the phenomenon of an IR 
saturation mechanism in which the extreme IR strength of the many contributing quark-quark 
couplings is easily modelled by this one effective gluon correlator. This was used in the context 
of the truncated Dyson-Schwinger (tDSE) approach by Munczek and Nemirovsky jllj who used 
a delta-function form. Of particular dynamical interest is the comparison of the GCM D liv {p) 
with one constructed theoretically from only a gluon correlator and vertex functions, say from 
continuum or lattice modellings, for this gives some insight into the IR strength of the higher order 
gluonic couplings. 



3 Hadronisation 

Hadronisation is a generalisation of the bosonisation concept to include the fcrmionic (baryonic) 
states. Bosonisation is naturally and conveniently induced using the FIC methods, and indeed the 
GCM is a spinor-boson field theory which may be exactly bosonised in 4-dimensional spacetime. 
A key feature of the GCM bosonisation is the use of bilocal fields which simply indicates that the 
complete theory is expressible using two-point correlators as an equivalent set of functional integra- 
tion variables. As we shall see the bilocal bosonisation (which precedes the hadronisation) not only 
generates the mesonic effective action, but also carries structural information which is essential in 
understanding the couplings of these mesons. This information (i.e. vertex functions) is essential 
to the prediction of hadronic properties and interactions. If we throw away this emergent infor- 
mation we are left with a non-renormalisable effective field theory with no predictive properties. 
There are two alternative first stage bosonisations of the GCM. The first in 1985 in which Cahill 
and Roberts|jJ introduced colour singlet and colour octet mesonic correlations had two problems: 
(a) the colour octet correlations did not appear to have any clear physical significance, and (b) the 
baryonic states were not manifest. However in 1989 Cahill, Praschifka and Burdenjl2| discovered 
a meson-diquark bosonisation which was induced by new colour and Dirac-spinor Fierz identities 
and which involved the colour singlet mesonic states and the colour 3 (qq) and 3 (qq) diquark/anti- 
diquark states, these being the very colour subcorrelations present in colour singlet baryons and 
anti-baryons. The meson-diquark bosonisation then led to the meson-baryon hadronisation of the 
GCM (Cahill 1989) |l3|]. This hadronisation automatically produced the covariant Faddeev for- 
mulation of the constituent baryon states, and also the meson dressing of these states; see later 
sections. In this review we concentrate mainly on this modern hadronisation of the GCM, which 
is summarised by the following sequence of FIC transformations: 

Z = J DqDqDAexp(-SQCD[A,q,q\+r}q + qr}) 

« J DqDqDA exp(-SGCM [A, q, q]+rjq + Vl) (GCM) 

= J DBDVDVcxp(-S[B,V,V}) (bilocal fields) (9) 
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= J DNDN..DnDpDuj...exp(-S h ad[N,N,..TT,p,u!....]) (local fields). (10) 

The basic insights are that (i) the quark-gluon dynamics in (0) is fluctuation dominated, 
whereas the hadronic functional integrations in (10) are not, (ii) the bilocal stage in (0) produces 



the constituent quark effect as the dominant configuration, and (iii) this entails the IR saturation 
effect and the dynamical breaking of chiral symmetry and its significant consequences, and (iv) 
the induced hadronic effective action in ( [Tof ) is nonlocal. The hadronisation has also been further 
studied in @ [0| [l| ([§ |l| 

3.1 Meson-Diquark Bosonisation 

Here we review the meson-diquark bosonisation 0j of the GCM, giving (||). The FIC techniques 
amount to analogues of various 'tricks' of ordinary integral calculus. Integrating out the gluon 
fields Z becomes, 

Z[V,V]=[ D QDqexp(-S[q,q}+r}q + qr^j, (11) 
where 

S[q,q] = J(q{x){ 1 .d x + M)6\x-y)q{y) + \q{x)^q{x)Df v ^ (12) 
Using the new Fierz identities [12] the quartic term in (|l2|) is rearranged to give 



S[q,q] = d 4 xd 4 y q(x)j.d6 4 {x - y)q(y) - -q{x)^-q{y)D(x - y). 

«(*) - -q{x)^q{yf T D{x - y)q{y)<* ^±q{ X )\ , (13) 

with q c — Cq , q c = qC. The Fierz identities are the two Dirac matrix identities 7^7*1 = K^ u Kf s 
where {K a } = {l,i l5 , -fc-f, ^7 P 7 5 }, and 7 £ s7 £ 1 = (K a C T ) rt (C T K a ) us where C = 7 2 7 4 , C 2 = 

-1 and Cy<C = 7^ T . We also use the set {K*} = {1, -z 75 , ^7 M , 2727^75}, then irfjfiT 6 ] = 
■iS ab . For the colour algebra Q A^A^ = ftW^-y + § I]p=i e pa7 e P<5/3> wnne for the iV/ = 3 
flavour algebra, SijSki = F^F^- for the mesons where {F c ,c — 0, ..8} = {-73I, 775, ■■, 4g} and 
= ff/ fe ff/. for the diquarks, where {H* , / = 1, ..9} = {F c , c = 7, 5, 2, 0, 1, 3, 4, 6, 8} and where 
are the generators of SU(S) in the usual Gell-Mann representation. We define the tensor 
products {M^} = W^K a F c } and {M$} = {i^/lK a ePH^}, where (e p ) Q/3 = e paf3 . We see that 
q(y)M^ n q(x) are l c bilocal fields with the flavour (1/ or 8/) determined by the flavour generators 
({F } or {F 1 '-" 8 }) in M^, while q(y) cT M$q(x) are 3 C bilocal qq fields with the flavour (3/ or 6/) 
determined by the flavour generators ({H 1,2 ' 3 } or {7J 4 '- ' 9 }) in Aft. These results follow from the 
colour and flavour representation of the quark fields. The (integral) spin of these boson fields is 
determined by the K a . By rewriting (O) as (|l^) we can initiate a bosonisation which is adapted 
to the attractive channels implicit in (]lj). The 1985 GCM colour Fierz identity Q leads to the 
colour 8 channels which are repulsive. 

We make the first FIC change of variables by noting that the quartic terms in exp(— S) may 
be generated by the following bilocal FIC integrations, 



J DqDqDBDVDV* exp ^ J -q(x)(j.d + M)5 i (x - y)q(y) 



B e (x,y)B e (y,x) V^x,y)V^(x,y)* 
2D(x - y) 2D(x - y) 
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]y[0 i yr4> 

-q{x)^q{y)B e {x,y) - -q( x )^±q(y)^ {x,yf 



V*(x, y)q{yf T ^±q(x)] + I \rjq + q V )) , (14) 



1^, - , ,^Ml 

where B (x, y) = B e (y, x)* are 'hermitean' bilocal fields. Integration over the quark fields completes 
the change of variables to bilocal meson and diquark fields, 



z iv,v] =J DBDT>DV*(DetJ : - 1 [B,T>,T5\)?.exp(J 



B 9 (x,y)B e (y, 



2D(x-y) 

f V+(x,y)V*(x,y)* 1 f^A 
J 2D(x-y) + 2j^ U )> 



(15) 



— —V r~ lT 

where 9 = (f},-r, T ), T~\B, V, V] = ( ^ ^ 



G-\x 1 y,B) = { 1 .d + M)5\x-y)+B{x 1 y), B{x,y) = B 6 (x,y)^, (16) 

V(x,y)=V*(x,yy^C T and V(x, y) = V+(y, x)C T ^. 
Using the determinant identity (l§] Detf- 1 = {Det^- 1 )) 2 Det(l + VG T VG), 
Z = J DBDVDV* cxp (rrLn{G[B]- 1 ) + ~TrLn(l + VG[B] T VG[B])+ 

QTe T \ . (17) 



2D J 2D 2 
3.2 Baryons 

The diquark sector of the meson-diquark bosonisation generates [Q the colour singlet baryon states 
of the GCM. Consider the diquark part of Z; 



J DVDV* exp (^TrLn(l + VG T VG) - J + J (J*V + V* J) 



where the bilocal diquark source terms facilitate the analysis, and in which the B dependence of 
G[B] will affect both the non-trivial dominant configuration and the mesons, and will provide the 
meson-baryon couplings. The expansion 

TrLn{l + VG T VG) = 2J Tr{VG T VG) n , (18) 

n=l U 

gives single loop processes (Fig.||(a)) with the quark lines alternating in direction, in accord with 
their coupling to the diquark and anti-diquark fields. Using (|l^) the diquark part of the action 
has the expansion S[V*,V] = S n [V* ,V] and with S x = J V^{lSr d l )^V^ and S' = S - S u 

Z = exp(-S'[^, J DVDV*exp(- Jv*A- 1 V + J(rV + V\J) 

= eMS'i^, A]) exp (-TrLn(A^) + J J*A d j) . (19) 
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(a) (b) 

Figure 4: (a) Example of diagrams from expansion of diquark TrLn in ([L8|). (b) Diagrams after diquark 
field integrations. 

Keeping only the translation invariant part of B, BcQ(x-y) - later to be identified as the constituent 
quark (CQ) effect, A^" 1 has eigenvalues and eigenvectors (diquark form factors) from 

<l>q (A-YHp,q;P)ri(q;P) = \k(P 2 )Ti(p;P) (20) 



(2tt)4 

and we have the orthonormality, completeness and spectral equations 
e J ^ ' k 

A*/( P ,q;P) = J2ri(p;P)\ k (PT 1 ri(q;Py. (21) 

fc 

Using completeness we construct the local-diquark-field F1C representation 

TrLn{A- 1 ) = ^ f d*x [ -gl/n(A fe (P 2 )) = £ TrLn{\ k {-d 2 )5\x - y)), (22) 
k J J ^ ' k 

exp(-Trin(A- 1 )) - / Dd k Dd$ exp(- £ / d k (x)* X k (-d 2 )S 4 (x - y)d k (y)) . (23) 

J k ^ 

Introducing local sources jk(X) = J d 4 Yd^xrf(x, X — Y)* J^{x, Y), so that 

but keeping only a single component of the scalar diquark to simplify notation, 

Z\j*,j] = exp(-5'[^, ^])cxp (-TrLniA- 1 ) + J f(X)X (-d 2 )- 1 j(X)\ • (24) 

Evaluating the effect of the functional operator we find that Z[0, 0] has the form exp(VF) where 
W is the sum of connected loop diagrams, with the vertices now joined by the diquark correlators 
Ao(P 2 ) -1 = (P 2 + m (-P 2 ) 2 )/ 2 , in which mo(P 2 ) is the diquark mass function, and with To(p; P) 
at the vertices. Of particular significance is the infinite subset of diagrams which will be seen to 
have the form of three-quark (i.e. baryon) loops (Fig.[|(a)). These come with a combinatoric factor 
of 2 (except for the order n=3 diagram) which cancels the \ coefficient of the TrLn. These 3-loops 
are planar for even order, but non-planar, with one twist, for odd order. To exhibit the three-quark 
loop structure we show, in Fig.||(a) , a typical diagram from Fig^(b) after deformation, revealing a 
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closed double helix: a diagram of order n is drawn on a Mobius strip of n — 1 twists. This infinite 
series may be summed as the diagrams are generated by the kernel K, defined by the one- twist 
diagram, shown in Fig.||(b). 



(a) (b) 

Figure 5: (a) Similar to Fig^(b) after redrawing to reveal baryon loop. These loop functionals determine 
the core baryon mass spectrum, (b) Kernel of the nucleon loop. 

The weightings are such that all the double helix diagrams may be summed to TrLn(l + K) — 
TrK = Wb — TrK. Thus Z[0, 0] = exp(W^ + D R ), in which D R is the sum of the remaining 
connected diagrams. To determine the content of Wb we consider the eigenvalue problem (1 + 
K)^ — which, for B — Bcq, has the following momentum space form, and is illustrated in 
Fig| 

/ ^ K ^ p t 3 i P >t^ p ) = W p2 ) - V*%rton (25) 

K(p, q; P)M$ = Yl ^oe 7a seif l ij 5 C T G T C T t l5 ef 3Sp e Jlh T G\a 1 . (26) 

iS 





Figure 6: The core baryon eigenvector equation. The diquark vertex is shown shaded. The diquark is 
essentially the same size as the baryon. 

Equation (^BJ) is a bound state equation for a three-quark state in which the paired quarks 
form a scalar diquark. ^ in ( ^5| ) is a Dirac spinor and, as well, a rank-2 tensor in both colour and 
flavour, ft should also be clear why the baryon loop functional consists of twisted quark-diquark 
lines; (^) is precisely of the form of the Faddeev equations of standard three-particle theory - 
the twisting merely arose so that the diagrammatic representation of (^5|) in FigJ| has the same 
legs in the same positions on both sides of the equation; conventionally Fig.^ would be drawn in 
untwisted form. In fact ( pff ) was also derived by conventional non-FIC three-particle methods [Q] . 
See [|o[ |2lJ] for early quark-diquark models. We note that because (|2l| ) is a discrete sum, (^) is of 
separable form. We separate colour and flavour multiplets in j25| ) , by decomposing according 
to (3 <g> 3) c ® (3 <g) 3) F = (lc ® If) © (lc <8> 8 F ) ® (8 C ® If) (S c ® 8 F ), 
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+ l^t - \nl^ih + - ln l h s P1 - \*pi h + l*AAh]- (27) 

Each member^ = [....]) of one multiplet is then seen to be an eigenvector of 

/ Wr [K(p ' 95 p) * (<?; p) = (A(p2) - m{p] n (28) 

N\m] ,P p P . ,P q P 

K(p, q; P) = ^r ( - + q + |; - - p)^ +p + l---q), 

M(j-q) 2 r 1 G(-q-p)G(^+ q ), (29) 

where the N[m] depend on the multiplet, and we find N[lc <8> If] = —2, N[lc <£> 8f] = +1, 
N[8 C ® If] = +1 and N[8c ® 8f] = — |. The l c ® l f and 8 C ® 8f multiplets have negative 
values for iV and thus the quark rearrangement process is repulsive, and the corresponding A(— M 2 )s 
have no zeros. However the colour octet - flavour singlet 'baryons' have an N value which means 
they are degenerate in mass with the colour-singlet flavour-octet baryons. Like the diquark states 
we expect colour non-singlets to not have a mass-shell, i.e. to be confined. It is known that 
crossed gluon processes can perform this task for the diquarks Q (see sect.^), but there has been 
no corresponding analysis for these coloured baryonic states. 

Let us now construct, for the colour singlet states, an appropriate FIC representation for 
exp(TrLn(l + K)). To this end note that an eigenvalue for positive energy solutions, with degen- 
eracy 2 (for spin \ and |), has the form A^(F 2 ) = (M(P 2 ) + i^/P 2 a{P 2 ))F , (define F so that 
a = 1 when A = 0, then M k (P 2 ) are baryon running masses), while for negative energy solutions 
(anti-baryons) Al^ = (A+)*. Thus, from the spectral representation for 1 + K, 



exp(TrLn(l + K)) = exp (^Tn J dAx J-^j 



jMA^n + MAlV 2 ) 2 )] (30) 



where k sums the ground state and excited baryons states of ( |25| ) , the squares in the In terms arise 
from the spin degeneracy and n = 8 arises from the flavour degeneracy (the other baryon states 
are not shown). Therefore, with a — 1 for simplicity, 



exp(Trin(l + K)) = exp 11 J ' <^ x J 
= exp ( TrL 

V k 



d 4 P 



(P 2 + M k (P 2 f) 2 F^ 



( r d + M k (-d 2 ))F^(x-y) 



k 



DN k DN k exp - V / d 4 xiV fc (x)( 7 .9 + M k {-d 2 ))F 2 N k {x) , (31) 



in terms of N k and N k , each of which is a flavour octet of local baryonic spin | FIC variables. 
Hence the exponentiated sum of the closed double helix diagrams is representable as a (free) 
baryon field theory. The F k may be absorbed with a re-definition of the baryon fields. Other 
more complicated (including baryon multi-loops) diagrams are present and constitute a wealth of 
dressings and interactions between these (bare) baryons. 

3.3 Mesons 

We now briefly indicate how the non-diquark part of fll7|), S[B], gives the meson sector. The 
complete Z has the form 

Z= f DBexp(-S[B]- ]T TrLn(X k (-d 2 ; [B CQ ])S\x - y)) 

diquarks 
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+ ^^{^.8 + M k {-d 2 -[Bc Q ]))5\x~y)) + ....J. (32) 

baryons 

We first determine the dominant configuration (and constituent quark effect) Bqq 1 as the solution 
of the Euler-Lagrange equations S jgw^ = 0, which gives 



B 6 CQ (x,y) =D(x-y) 



tr(G(x,y,B CQ )^) 



(33) 



a non-linear equation for the {B 6 C q\. - a Dyson-Schwinger type equation, where '+....' are the 
diquark and baryon parts. Dynamically it describes the extensive self-energy of the quarks due to 
dressing by gluons, leading to the quark running mass. This is finite and extends over distances 
comparable to hadronic sizes and with an energy density that implies that it has the dominant role 
in determining hadron structure. Only colour-singlet translation-invariant solutions (depending 
only on x — y) are known. Fourier and inverse Fierz transforming ([33]) we obtain, on retaining only 
the meson contributions, 



4 f 1 



B C q(p)=c- I j—r jDip-q)^ .^ _ | KA | g __^ 7M- ( 34 ) 



Here c = § from the meson-diquark bosonisation, while normal Feynman rules would give c = 1, 
as does the 1 — 8 GCM bosonisation [Q. That c < 1 in the 1 — 3 — 3 bosonisation indicates 
that some strength is generated by additional mechanisms in ( |33| ) involving meson and diquark 
processes, i.e. the 'book-keeping' is more subtle in this bosonisation, and is so to avoid double 
counting. At present we adopt the practice of using c = 1 until these additional processes can be 
included dynamically. This CQ equation has unique solutions when M. ^ and G has the form 
G(q) = [iA(q)q.^+M+B(q)]^ 1 . Expanding S[B] about its minimum gives S[B] = J2 n =o 2 3 S n [B], 
where S n is of order n in B and, for example, S2 = \ J B e (A^ 1 ) ^ B^ . Introducing bilocal source 
terms in (^21) we have, with S' — S — 5*2, and showing only the meson part, 



Z[J] = J DBcxp(-S'[B}- S 2 [B}+ J B e J e ) 

= exp(-5'[^]) J DBexp(- | ^(A" 1 )^ + J B 9 J e ), 
= exp(-S ,/ [^-]) J Dm k exp(~y^ ^ J m k (x)\ k (-d 2 )m k (x) 



jkm k ). 



Here we have used techniques similar to that for the diquarks and {m k (x)} is an infinite set of local 
meson fields. Each m k corresponds to one physical meson type, and the are the eigenvalues 
of the meson form of (^)- a Bethe-Salpeter equation, which also gives the meson form factors 
T k (p 7 P). Applying the functional operator exp(— S'[jj]), 



Z = J Dm k cxp (^Z^2 \ J m k{x)\(-d 2 )m k {x) ~ S'[m k ]j 



(35) 



By evaluation of S'[m k ], and identifying the mesons by their quantum numbers, we obtain the full 
local FIC representation for the meson sector of QCD 23, p4| , as summarised in sect. |3.5| . 



3.4 Hidden Chiral Symmetry 

When the quark current masses M. — > the fundamental action S\q, g,^] has the important 
additional global Ul(Nf) ® Ur(Nf) chiral symmetry. The consequences of this follow naturally 
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through the GCM hadronisation. The first significant result is that the dominant configuration 
(HJ) has degenerate solutions |g| g|[ 

G(q; V) = [iA(q)q. 1 + VBiq)}- 1 = ?G(q; (36) 

where £ = y/V, V = exp(i v / 2757r a F a ) and {tt "} are arbitrary real constants |7r|<E [0, 2tt]. Thus in 
the chiral limit the dominant configuration is degenerate and is the manifold (Ul^Ur)/ H (a coset 
space) where H — Uy C Ul <£> Ur. Thus the chiral symmetry is represented as a hidden symmetry. 
This occurs because the action in (^2|) has the form of a 'Mexican-hat' in the relevant variables H. 
The Nambu-Goldstone (NG) bosons form homogeneous Ricmann coordinates for this dominant 
configuration manifold. There is a technical complexity in expanding S"[mfc] in the chiral limit 
which is caused by the action having degenerate minima, since we do not have a unique minimum 
about which to expand. First we need fields adapted to the compact dominant configuration 
manifold. For this we use the angles {it} as new field variables {tt(x)} |25|, |2(| in place of some of 
the B e (x, y). The dependence of the action on these dominant configuration variables is such that 
the action will increase only if the dominant configuration point is different at different space-time 
points, and so a derivative expansion in d^V{x) must arise. The Dirac algebra allows finally the use 
of the matrix U{x) = exp(iy/2TT a {x)F a ) where V{x) = P L U(x)^ + PrU(x) = exp{iy/2j 5 ir a {x)F a ). 
Then the NG sector of S[B] is 

f 2 

(fx^trid^Ud^) + Kl tr(d 2 Ud 2 U^) + n 2 tr([d^Ud^} 2 ) 

+ KsbridJJdjfldJJdjfl) + |tr([l - ^L]M) + .... (37) 

where f n) Ki,.. are given by explicit integrals [^5], in terms of A(q) and B(q) and p is the 
quark condensate parameter. In the chiral limit it is important to note (5[ |25|] that B(q) in the 
quark correlator is also the NG boson on-mass-shell form factor, T 7r (p;P = 0) = B(p). In the 
chiral limit the ground state pseudoscalars play a dual role: they are at the same time both the 
NG bosons associated with the hidden chiral symmetry and also qq bound states. Under a chiral 
transformation we find H C7(a;) -> U L U(x)U R . This is a derived result of the FIC analysis which 
is usually assumed in phenomenological modelling. We have included the lowest order term which 
depends on A4, i.e. for small breaking of the chiral symmetry by the quark current masses. 

The coupling of the baryon states to the above mesons requires us to keep the full B in 
analysing the baryon sector, and not just the dominant configuration value Bcq- However the 
long wavelength limit of the NG-boson-baryon coupling may be inferred from the chiral invariance 
of (02|). Now 



TrLn 



{ 1 .d + M{-d 2 ))8 A {x-y) 



— TrLn 



{ 1 .d + VM{-d 2 ))5 4 {x^y) 



reflects that invariance in (|32|), where V = exp(« v / 2757r a T a ), with {T a } the generators of SU(3f) 
8 representation. 

3.5 Hadronic Laws 

Gathering the above results and keeping only the low orders in the hadronic variables and in the 
derivatives (appropriate to a low-energy long-wavelength expansion) 



Z = J D7rDpDuj..DNDN...exp(-S had [n,p,Lj,..N,N..}) , 

Shad [ir,p,u,...,N,N,..] = 

J d 4 xtr{N( 1 .d + M + AM - M^2i^ a T a + ..)N} 



(38) 
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^[(^tt) 2 + mln 2 ) + J -£-[-p^- d 2 )p„ + (d,p,) 2 + mX] 



P 

-^-[p — ► w] - fpf^gp-K-Kp^ x d^-K - if u f„€n V<TT u) li d v n.d <T ir x d T ir 

^ feu fp fn Gijj Cfj, v a T UJ^ 0^ p a . 9 T 7T 

Xi 



80tt 



tri-K.FdpTT.FdvTr.Fdvir.Fdrir.F) 



2 c /^^crr vi V"-- 1 - w fi 



(39) 



in which the baryon octet is finally written as a rank-2 tensor, N = N a T a , where the {T a } 
are generators of the SU(3f) 3 representation. We have written Xj(P 2 ) = (P 2 + mj(P 2 ) 2 )f 2 
where rrij(P 2 ) are the running meson masses, but only the physical masses (from X(P 2 ) = 0) 
are shown above. The imaginary terms in this meson action are the chiral anomalies of QCD, 
including in particular the Wess-Zumino term. In (|3^) we have shown m v and AMq which are 
mass terms from the chiral symmetry breaking quark current masses, while Mq is the 'chiral mass' 
of the constituent baryons (see sect.||). For non-zero quark current masses the NG boson masses 
{771^} and the baryon octet mass splittings {Amo} are seen to satisfy the Gell-Mann-Okubo and 
Coleman- Glashow formulae p7|| . 

In general the coupling terms in the hadronic action are non-local and the actions in ( |37| ) and 
( |39| ) will also contain higher order derivative terms like tr{N (mjv \f2i"fr > d 2 ir a T a + ..)N}. These 
should not be thought of as 'different' couplings, but rather as just arising from the expansion of 
the meson-baryon vertex function Tq(p, q). Hence rather than making this effective action non- 
renormalisable, as often occurs in effective actions, such terms when properly retained as parts 
of complete vertex functions actually render loop diagrams finite (an example is the pion-nucleon 
loop in sect||). 

The full non-local meson sector of (|3^) |p3 has been used in many studies, such as the p — * mr 
decay p3| , u> — p mass splitting p8[ , charge symmetry breaking via p — lo mixing j2j| , pion loop 
contribution to p — lo mixing ]3C| and pion and p meson observables |3lJ , for extensions to include 
electromagnetic interactions ]32| , |33| ] , and for rj and rj' p(| . For the pion loop contribution to the 
electromagnetic pion charge radius see |34|] . The chirally invariant form of the NG boson sector in 
(p7|) in particular has been investigated in pBl and the it — it scattering lengths in p6| . 



4 Effective Gluon Correlator 

We now consider the detailed implementation of the GCM. First we must determine the effective 
gluon correlator by fitting GCM observables to experimental data. This involves the determination 
of the dominant configuration, i.e. the constituent quark effect. This effective gluon correlator and 
quark correlator are then used in the constituent-meson BSE equations in order to determine meson 
masses, and also f^. The dominant configuration is defined by equations, 

* S 1 ,=0. (40) 



SB(x,y) 



Of the set Bcq{x,d) only A(x — y) and B(x — y) (their Fourier transforms appear in (42)) are non- 
zero bilocal fields characterising the dominant configuration. They are also translation-invariant. 
This is the dynamical breaking of chiral symmetry. Such non-zero dominant configurations are 
also known as condensates. Writing out the translation invariant CQ equations we find that the 
dominant configuration is indeed simply the constituent quark effect as they may be written in the 
form of @ or (||), or 

G- 1 ( P )=ik + ™+^J ^D MV (p-q)<y ll Giq)<y v , (41) 
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which is the gluon dressing of a constituent quark. Its solution has the structure 

G(q) = (iA(q 2 ;m)q.j + B(q 2 ;m) + m)~ l = -iq.-f<r v (q 2 ;m) + <r s (q 2 ;m). (42) 

In the chiral limit there are more Bcq fields that are non-zero, and the resultant degeneracy of 
the dominant configuration is responsible for the masslessness of the pion. The constituent quark 
correlator G should not be confused with the complete quark correlator Q from (||) which would 
be needed to analyse the existence or otherwise of free quarks. The G on the other hand relates 
exclusively to the internal structure of hadrons, and to the fact that this structure appears to be 
dominated by the constituent quark effect. The evaluation of Q is a very difficult task, even within 
the GCM, while G is reasonably easy to study using ([4l|). The truncation of the DSE in which 
the full quark Q is approximated by this G amounts to using a mean field approximation (see 
(|53|)); however from the tDSE there is no systematic formalism for going beyond the mean field as 
there is in the GCM. The hadronic effective action in (^) arises when S[B, ..] is expanded about 
the dominant CQ configuration; the first derivative is zero by ([lO]), and the second derivatives, or 
curvatures, give the constituent or core meson correlators G m (q,p; P) 

G ^ P ^ = F - T \ sB(J)!b(u,v) I**)' ™ 

after exploiting translation invariance and Fourier transforming. Higher order derivatives lead to 
couplings between the meson cores. The G m (q,p; P) are given by ladder-type correlator equations 



and the non-ladder effects are inserted by the final functional integrals in (10), giving the complete 
GCM meson correlators Gm(q,P', P)- It is interesting to note that the truncated and modified DSE 
in Maris and Roberts |37) are identical to (|4l| ) and (fl3j), in the form of (HJ). 

In the fitting to meson data the u> and ai mesons are described by these constituent meson 
correlators; that is, we ignore meson dressings of these mesons. The mass M of these states is 
determined by finding the pole position of G m (q,p; P) in the meson momentum P 2 — —M 2 and 
this, or equivalently the meson version of (20) which for the mass-shell has A(P 2 ) = 0, leads to the 



homogeneous vertex equation 

l>; P) = ~\f (Si^Ws - PhMi + f )r(g; P)G(q - £) 7 „. (44) 

To solve pl| ) for various D fJil/ (p) and then to proceed to use A(s) and B(s) in meson corre- 
lator equations for fitting observables to meson data is particularly difficult. A robust numerical 
technique is to use a separable expansion^) as follows. We have in the Landau gauge 

D^(p) = (6^-^)D(p 2 ), and ^(p) = - ^)G(p 2 ), (45) 
pr pr 



where D(p ) = g{p )G{p )g{p )■ First expand D{p—q) in (41) into 0(4) hyperspherical harmonics 



D(p-q) = D Q (p 2 ,q 2 ) + q.pD^p 2 ^ 2 ) + (46) 

Do(p 2 ,q 2 ) = - / dfisin 2 (3D{p 2 + q 2 - 2 M cos/3), ... (47) 
I" Jo 

Introduce a multi-rank separable D expansion (here n = 3) 

D (p 2 ,q 2 )= Yl r *(P 2 ) r ^ 2 )> ( 48 ) 

i=l,n 

and the constituent quark equations then have solutions of the form (in the chiral limit) 

B(s) = ]T Bi(s), B^s) = biTiia), a s (s) = £ a a (s)i, (49) 
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b 2 = 4n 2 f sdsBiisJasis), and B^a) = . §^ . ,„ ■ (50) 

However rather than specifying r, in we proceed by parametrising forms for the a S i and 
a v ; the then follow from (pi) and (S3): 



cr s (s)i = ciexp(-dis), ct s (s) 2 = c 2 . 



2s 2 -rf 2 (l-exp(-2g 2 /rf 2 )) x 2 
2^4 



/ \ ^ 2/( S )-d 3 (l-exp(-2/( S )/d 3 )) V 2 1/2 

ff s (s)3 = c 3 I I , /(s) = s{ln(T + s/A )) ' 

2 S -/? 2 (l-exp(-2 S //3 2 )) 

<M«) = ^2 • ( 51 ) 

The three a S i terms mainly determine the IR, midrange and UV regions; the <r s (s)3 term de- 
scribes the asymptotic form of a s (s) ~ 1/s 2 ln(s/A 2 ) for s — > oo and ensures the form for 
D(s) ~ l/sln(s/A 2 ). The parameter r = 10 ensures that 0-3 is well behaved at small s, but 
otherwise has no effect on the fit. With these parametrised forms we can numerically relate, in 
a robust and stable manner, the parameter set in Table 1 to the mass of the ai(1230MeV) and 
w(783MeV) mesons from ([5|), to 7^(93. 3MeV) and to experimental a(s) points (see insert in Fig.Q) 
from the Particle Properties Data Booklet for s > 3GeV 2 . 

The translation invariant form for the effective gluon correlator is easily reconstructed by using 
D{p 2 ) = Dob 2 ,0) in (Eel) and then from @ 



[P > ^ b\ a s (0) 2 p 2 a v (p 2 ) 2 + a s (p 2 f ' [ ° 1 

With the parameter set in Table 1 the resulting quark-quark coupling correlator D(p 2 ) is 
shown in Fig.0 and Figj8[ A significant feature of QCD is that the infrared dominance, as revealed 
by the large value of D(s) at small s, causes the CQ equations to saturate, i.e. the forms of the 
solutions A(s) and B(s) at low s are independent of the detailed IR form of D(s). This saturation 
effect means that low energy QCD is surprisingly easy to model, and this effect is utilised in the 
GCM. 

In Fig.||we plot the new GCM98 quark-quark coupling correlator D(s) which shows excellent 



overall agreement with the Jain-Munczek D(s)|38, p9|, and we also plot D(s) = g(s)G(s)g(s) 



constructed from the Marenzoni et al. El] Q{s) and Skullerud p2[ g(s) lattice results which shows 
agreement down to s — 1.8GeV 2 . The normalisation and shape of the Marenzoni et al. Q(s) 
agrees with that of Leinweber et al. pC[ | . However the normalisation of the much more difficult 
lattice computation of g(s) is uncertain and we have chosen it so that the combined lattice D(s) 
agrees with the experimental Particle Properties Data Booklet for s > 3GeV 2 . As shown in Fig.|| 
all three D(s) depart from the two-loop form below s — 2.5GeV 2 . The difference between the 
GCM98 (or Jain-Munczek) and the lattice construction could be an indication of contributions to 
the quark-quark coupling from higher order gluon self-couplings at low energy, since processes like 
Fig.y](c) would be included in the GCM fit, but are not in the lattice construction. The earlier 
GCM95 and GCM97 gave D(s) that differed mainly in the asymptotic region, see plots in |To|| . 

The insert in FigJ| shows the g(s) from g 2 (s) = D{s)/Q{s) that then follows from our analysis. 
This is the effective quark-gluon coupling vertex if the gluon correlator is taken to be that of 
Marenzoni et al. (or Leinweber et al.) and here the error bars now indicate combined errors and 
uncertainties from the lattice spacing. Also shown is g(s) from Skullerud |f2|| with the normalisation 
as discussed above. 

We summarise in Table 2 some of the hadronic observables that may be computed within 
the GCM, showing in particular their sensitivity to the evolving modelling of the effective gluon 
correlator. Further meson observables are reviewed in Tandy 0. 
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Figure 7: Plots of D(s) (GeV~ 2 ). Solid line is GCM98; shortdash line is GCM95; longdash line is GCM97. 
Data plot is lattice pure-gluon Q{s) from Marenzoni el at, and so has no quark-gluon vertex. Insert is fit 
of GCM98 to a(s) of Particle Data Book. 



Table 1: GCM1998 <r s (s) and a v (s) Parameters. 



Cl 




di 


3.620GcV- ;i 


a 


0.4956GeV 


C2 


0.0281GeV 7 




1.516GcV 4 


A 


0.234GcV 


C3 


0.0565GeV 3 


d 3 


0.7911GeV 2 
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s (GeV 2 ) 

Figure 8: Plots of D(s) (GeV -2 ). GCM98 is the upper solid line which is almost identical to the Jain and 
Munczek D(s) (dashed line); lower solid line is two-loop form with A = 0.234GeV, Nf = 3 and r = 10; 
data plot is combined lattice data for g(s)G(s)g(s) with G(s) from Marenzoni et al. (as in (a)) and g(s) 
from Skullerud. Insert shows g(s) from Skullerud (lower data plot), and from GCM98/(J(s) (upper data 
plot). 
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Table 2: Hadronic Observables. 



\ 1 \~\ coi'^/n \~\ 1 o 

VJUBcl VdlLIlt; 


V 1 V - 1VX 1 :>:>•} 


PfMI QQ7 

V l V . \ 1 1 .7.7 l 


V l V . \ 1 1 . 7 . 7 ( ~ 


1 . A | H / 1 1 1 ( ( J 1 V 


f 


QQ OA/In A/* 

yo.uivie v 


QQ OMnV* 

yo.zivie v 


Q9 zLni\/ToV* 
yz.Ttuivie v 


QQ Q A/To A/ 

yo.oivie v 


d\ meson itlrss 


izouivie v 


1001 MpV* 

izonvie v 


1 93QA/IW* 

izoyivie v 


1 93flA/TnV 
izouivie v 


7T meson mass(for xxiu.d) 


loo.oivie v 


1 38 ^A/ToA/* 

loo.oivie v 


1 38 ^.A/ToV* 
loo.oivie v 


1 38 ^A/ToA/ 
loo.oivie v 


a(s) 






lllieLl 


T 


TC lnp'jnn in 3 kci ( for tn 1 


496MeV* 






496MeV 


(m u + m d )/2\ R (u = lGeV) 

\ Li ' Lb / J | il, \ i J 


6.5MeV 


4.8MeV 


7.7MeV 


«8.0MeV 


rrie\n(u — lGeV) 


135MeV 






130MeV 


uj meson mass 


804MeV 


783MeV* 


783MeV* 


782MeV 


Oq 7r — 7r scatt. length 


0.1634 


0.1622 


0.1657 


0.26 ± 0.05 


Oq 7r — 7r scatt. length 


-0.0466 


-0.0463 


-0.0465 


-0.028 ± 0.012 


a\ 7r — ir scatt. length 


0.0358 


0.0355 


0.0357 


0.038 ± 0.002 


a§ 7r — 7r scatt. length 


0.0017 


0.0016 


0.0018 


0.0017 ± 0.003 


a\ it — 7r scatt. length 


-0.0005 


-0.0005 


-0.0003 


.00013±0.0003 


r,r pion charge radius 


0.55fm 


0.53fm 


0.53fm 


0.66fm 


i (O^)nucleon-core mass** 


1390MeV 


1435MeV 


1450MeV 


>1300MeVff 


const, quark rms size 


0.51fm 


0.39fm 


0.58fm 




chiral quark const, mass 


270MeV 


267MeV 


325MeV 




+ diquark rms size 


0.55fm 


0.55fm 


0.59fm 




+ diquark const, mass 


692MeV 


698MeV 


673MeV 


>400MeV 


1+ diquark const, mass 


1022MeV 


903MeV 


933MeV 




0~ diquark const, mass 


1079MeV 


1049MeV 


1072MeV 




l - diquark const, mass 


1369MeV 


1340MeV 


1373MeV 




MIT bag constant 


(154MeV) 4 


(145MeV) 4 


(166MeV) 4 


(146MeV) 4 


MIT N-core (no cms corr.) 


1500MeV 


1420MeV 


1625McV 


>1300MeVft 



* fitted observable; - not computed or not known; f a(s) from Particle Properties 
Data Booklet; GCM1995: §; GCM1997: §; GCM1998: @. 



** only + diquark correlation; 1 + diquark correlation lowers nucleon core mass, 
ff nucleon core mass (i.e. no meson dressing). 



5 Constituent Quarks 

The constituent quark effect jl3| is the dominant effect in determining the structure of hadrons, 
and also their response in scattering events, particularly deep inelastic scattering. The constituent 
mass effect manifested itself in the early studies of baryon magnetic moments. The GCM analysis 
reveals rather directly both the effective mass and the effective size of the constituent quarks, 
and relates these to the effective gluon correlator. We consider chiral limit constituent quarks 
and we carefully define constituent quarks and constituent hadrons as those constructs appearing 
in exponentiated effective actions in the functional formulation of the GCM and distinguish them 
from the exact correlations, which follow from the complete functional integrations. We can express 



the full quark correlator Q(x,y) in terms of the bosonised FIC variables, with (J2J) and (17) giving 



_ J DBDVDV*G(x, y, g)exp(-S[g, V, V*}) 
y[X,V) J DBDVDV*cxp(-S[B,V,V*}) ' ( ' 

where G{x, y, B) is defined in (fl6|). The bilocal field functional integrals can be further decomposed 
into local hadronic functional integrals. The constituent quark effect appears as the dominant 
configuration about which the meson-diquark bosonised GCM action in ( p3[) (and elsewhere) is 
expanded. This is characterised by T> — and two of the B ^ (the Bcq)- At the dominant 
configuration G{x,y 1 Bco) is the constituent quark correlator which appears in the constituent 
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meson and baryon correlators. The Bqq eqns.((34j) have the form (we use a Feynman-like gauge to 
simplify the discussion, but the more realistic Landau gauge may be used), 

^> = f/(^-« WprW (54) 

ri/ „ n , 8 f d A q , s A(q 2 ) , s 

[A{p a )-^ = 3/ ^ w ^-*?WtWf- < 55 > 

Here B(q 2 ) and C(q 2 ) — A(q 2 ) — 1 are the Fourier transforms of the two BcQ(x,y) fields, and 
where translational invariance is used: Bcq(x,u) Bcq(x ~ y)- We see from (|53|) that the full 
quark correlator is given by the dressing of the constituent quark correlator by various hadronic 
fluctuations. In the context of the tDSE (Q) and (^5|) (in Landau gauge) are often used as an 
approximation |}7| to the full quark correlator; Q(x,y) w G(x,u,Bcq), thereby confusing the full 
and constituent quark correlators. The GCM formulation clearly reveals that the tDSE approach 
is actually using a mean field approximation. The additional processes manifest in (|5^) show up in 
the FIC hadronisation as the dressing of constituent hadrons by other hadrons, of which the pion 
dressing of the nucleon is the most pronounced example. These dressings incorporate additional 
processes corresponding to further dressing of the constituent quarks as well as further interactions 
between the constituent quarks. 

When considering only the constituent quarks the bilocal effective action in ( p3[ ) may be 
simplified by keeping only the translation invariant B and C dependence. This leads to a reduced 
action per quark flavour, 



S CQ [B,C]=v(-^- 4 J dqqHn(A 2 (q 2 )q 2 + B(q 2 ) 2 ) 



-7r / dxx 3 v \ + — 7r / dxx 3 \ ; (56) 



4 J D(x) 2 J D{x) 

where V is the (infinite) spacetime volume. The minimization of Scq[B,C] gives ( |54| ) and (p5|). 
The action Scq[B, C] has the form of a sum of a kinetic energy term (defined as that part which 
is local in momentum) and a potential energy term (defined as that part which is local in relative 



spacetime). Both have unconventional forms because (56) describes self-interaction effects. The 
kinetic energy term involves the constituent quark running mass M(s) — B(s)/A(s). 

Since the key constituent quark effective mass is associated with the kinetic part of (|56|), we 
subtract the B — form. In this way we compare the non-perturbative configuration with the 
perturbative configuration, and it is this difference which also generates the MIT bag constant 
discussed in sect .[101 



/ 12^ 2 /•+» f f A(q 2 ) 2 q 2 + B(q 2 ) 2 \] 9 2 /■+» \ 3 B(x) 2 ]\ 

The kinetic energy and potential energy integrands, indicated by the square brackets in (p^) , are 
shown in Fig^. Fig.^(a) also shows the quark running mass. The chiral limit constituent quark 
mass of approximately 300MeV is revealed as the value of the running mass that dominates the 
kinetic energy integrations. The width of the g-integrations being a 'fermi-motion' effect. The 
integrand of the potential energy term shows that gluon exchanges up to some 1.2fm are relevant. 
Hence we see directly that the constituent quark characteristics are implicit in the action (^6|). 
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Figure 9: (a) Quark running mass M(s) (solid line) and the integrand of the kinetic energy part of the 
constituent quark action, (b) Integrand of the potential energy part of the constituent quark action. 

6 Constituent Mesons and Diquarks 

From the hadronisation we saw that the constituent mesons and diquarks arise as ladder BSE 
states; with their correlators described as the generalised curvatures of the meson-diquark boson 
action when expanding about the dominant or constituent quark configuration. This provides a 
particularly instructive insight into the true role of such ladder BSE states. Traditionally these 
ladder states arose as a severe truncation of the full coupled DSE, however we see that they 
actually play a key dynamical role in the hadronisation in that they naturally arise as appropriate 
FIC variables, rather than from some unstructured approximation scheme. The crossed diagrams 
that are normally neglected in the tDSE are automatically inserted in the GCM via the hadronic 
functional integrations in (|38|), and it is these processes that convert the constituent mesons and 
baryons into the observable hadronic modes, as seen later in sectj^ for the nucleon. Hence the 
GCM hadronisation reveals a 'book-keeping' that was not previously known. In this connection 
Bender, Roberts and von-Smeckal have found evidence that these additional crossed gluon 
processes may be responsible for confining the diquarks, which is particularly interesting since the 
ladder BSE have a mass-shell for the diquarks, and it is these masses which are shown in Table 2, 
see also [Q. That study explicitly included the crossed diagrams by extending the BSE equation 
to include the lowest order crossed diagram. However in the GCM hadronisation such crossed 
diagrams are seen to arise from meson exchanges, as shown in Fig.^. Evidence from the nucleon 
computations suggests that due to the high 'constituent mass' of the diquarks the presence of the 
pole in the diquark correlator at the mass-shell is not dynamically significant. 
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(a) (b) (c) 

Figure 10: (a) Diagram shows a meson exchange from the functional integral in (pr^), dressing constituent 
ladder diquark correlations, (b) a low order gluon process within the meson exchange. In (c) we redraw 
(b) to show the crossed gluon processes inserted via meson exchange. 

The connection between the bilocal meson-diquark action and the ladder BSE states inspired 
an insightful alternative to the solving of the linear homogeneous BSE. This analysis in Cahill, 
Roberts and Praschifka (45| involved re-formulating the ladder BSE, in which the boson mass of 
interest M occurs implicitly in the quark correlators, as an explicit mass functional. As an example, 
for the scalar diquark state + we find 

Mfrl 2 --^i- [ A r(g)2 I 9 fd^^l (58) 

[L 1 /[rp J (2^)4 A( q )w + B{ q y + /[rf J d(x) ' [b *> 

where f[T] is a normalisation functional Ej|. The minimisation of M[T] with respect to the boson 
vertex function T then yields the mass M > 0, at least as a good approximation for the low mass 
mesons and diquarks. A similar expression for the pion mass functional M ff [r] is given by (|58|), 
but with 9 — > 9/2 in the second term (the Landau gauge version, in the near chiral limit, is given 
in (|73|)). Then the minimum of M,r[r] = = arises for T(q) = B(q), as follows from the 
minimisation of Scq [B] in ( |57|) . However another feature of the GCM comes into play when we 
use the property that all the low mass mesons and diquarks have approximately the same r(q), 
and that this is is merely the B(q) function of the constituent quark correlator (an exact result for 
the pion) . This intricate relationship is again a consequence of the dominance of the constituent 
quark effect. Hence using T = B in (|58|) then yields an explicit value for the constituent mass of 
the + diquark. Similar mass functionals for the other states are given in GjJ. 



7 Constituent Nambu-Goldstone Mesons 

The properties of the pion continue to be the subject of considerable theoretical and experimental 
interest in QCD studies. The pion is an (almost) massless NG boson and its properties are directly 
associated with dynamical chiral symmetry breaking and the underlying quark-gluon dynamics. 
The GCM is particularly effective in revealing the NG phenomena that follows from the dynamical 
breaking of chiral symmetry (2^, . Indeed the GCM results in the complete derivation of the 
Chiral Perurbation Theory (ChPT) phenomenology, but with the added feature that again the 
induced NG effective action is non-local, so that the usual non-renormalisability problems do not 



arise (see sect. 3.4). However again we must distinguish the full NG degrees of freedom from 
the constituent NG modes, as again this distinction is often missing, particularly in the tDSE 
formulation. The full NG (pion) correlator Q n is the connected part of 



Q v (x,y,z,w) =1 DqDqDAq(x)i'j 5 T i q(y)q(z)i'Y S T i q(w) exp(-SacM[A, q, q\) (59) 
or, from the hadronisation (fl0|), 

g T (X,Y) = [ DTr..DNDN...Tr(X)n(Y)cM-Shad[7r 1 ...,N,N,..]) (60) 
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in which X = and Y = are 'centre-of-mass' coordinates for the pion. We note that 



now the pion field appears in Shad[^, N, N, ..] in the exponent of (60), with an effective-action 



mass parameter m^; see (39). It is important to clearly distinguish this mass, together with the 
equations which define its value, from the pion mass that would emerge from the evaluation of the 
functional integrals in (^9|) or (^). Equation (^) or ( |60| ) defines the observable pion mass, whereas 
the mass in the exponent defines the constituent pion mass. There is no reason for these to be 
equal in magnitude, though they may well both be given by the generic Gell-Mann-Oakes-Renner 
(GMOR)|§ formula. 



m„ = 



1 -Ji\ (™«( s ) + m d {s))o s {q 2 ) + 0(m 2 ). (61) 

Recently there has been a derivation |f7], H [l9) of an alternative mass formula and the demon- 
stration of its equivalence to the GMOR formula for the constituent pion mass, which we briefly 
consider. That derivation exploited the intricate interplay between the constituent pion Bcthe- 
Salpter equation (BSE) and the non-linear CQ equation, resulting in the new expression 

"** = 7^ / ( £ -( s )^( s ) + se v (s)a v (s)) c(s) + 0(m 2 ), (62) 
where c(s) is the function 

where A(s) and B(s) are the chiral limit of A(s;m) and B(s;m), and where e s (s) and e v (s) are 
functions which specify the response of the constituent quark correlator to the turning on of the 
quark current mass, defined, for small to ^ 0, by the expansions in m(s) 

B(s; to) + m(s) = B(s) + m(s).e s (s) + 0(m 2 ), (64) 

A(s, to) = A(s) + m(s).e v (s) + 0(to 2 ). (65) 

For large space- like s we find that e s — > 1, but for small s we find that e s (s) can be significantly 
larger than 1. This is an infrared-region dynamical enhancement of the quark current mass by 
gluon dressing, and indicates the strong response of the chiral limit constituent quark correlator 
to the turning on of the current mass. A plot of e s (s) is shown in [^7). The GCM involves the 
solution of CQ integral equations ( pll| ) for the constituent correlation functions. Separating this 
CQ into its scalar and vector parts , we obtain in Landau gauge 

Til 2 \ a [ A n / ^ B(q 2 ;m) +m(s) 

B(p ;m) =4 / -—^D(p-q). — — (66 

J (27r) 4 q z A{q z ] m) z + {B{q z ; to) + to(s)) z 



[A{p 2 -m)-l]p 2 



d^q 



(27T) 



jD{p-q) \p.q + 



,<i-ip - q)p-{p - q) 



(p-qf 



A(q 2 ;m) 



q 2 A(q 2 ;m) 2 + (B(q 2 :m) + m(s)) 2 ' 



(67) 



which are ( p4| ) and (55) (which however are in Feynman gauge) with non-zero quark current mass. 
We have included, for generality, a phenomenological momentum dependent current mass m(s) for 
the quarks. The usual procedure is to introduce a cutoff A, and to choose an s independent to(s), 
but with the value of that constant to now A dependent. This is in fact our final choice. However 
our analysis supports a more general result with m(s) dependent on s; then the GMOR relation in 
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(|6l| ) has the m(s) included inside the integration. Using Fourier transforms ( p6|) may be written 
in the form, here for m = 0, 

4 <7 s (a:) 

which implies that knowledge of the chiral-limit constituent quark correlator determines the effec- 
tive gluon correlator. Multiplying ( |68| ) by B(x)/D(x), and using Parseval's identity for the RHS, 
we obtain the identity 

j4 B(x) 2 A f A 



^# a4 ]# B(W) " (69) 

The second basic equation is the ladder form BSE for the constituent pion mass-shell state, 
which arises from the mesonic fluctuations about the minimum determined by ( |66| ) and (|67]) , 

r'(p, P) = 2 J j^D(p - q)trs F (G + T°G-Tf)r°(q, P) (70) 

where G± = G(q±^). This is (Q) in Landau gauge for isovector NG bosons, and only the dominant 
r = T'TH^ amplitude is retained. The spin trace arises from projecting onto this dominant 
amplitude. Here {T b , 6 = 1, .., N 2 - 1} are the generators of SU(N F ), with tr{T*T<i) = §£ /9 . 

The BSE ( fr0| ) is an implicit equation for the mass shell P 2 = —m 2 . It has solutions only in the 
time- like region P 2 < 0. Fundamentally this is ensured by ( |66| ) and ( |67j ) being the specification of 
an absolute minimum of an effective action after a bosonisation. Nevertheless the loop momentum 
is kept in the space-like region q 2 > 0; this mixed metric device ensures that the quark and gluon 
correlators remain close to the real space-like region where they have been most thoroughly studied. 
Very little is known about these correlators in the time-like region q 2 < 0. 

The non-perturbative quark-gluon dynamics are expressed here in ( |66| ) and ([37]). Even when 
m = (|(36| ) can have non-perturbative solutions with B 7^ 0. This is the dynamical breaking of 
chiral symmetry. When m — (|7(]) has a solution for P 2 — 0; the Goldstone theorem effect. For 



the zero linear momentum state {P = 0, P4 = 0} it is easily seen that (70) reduces to ( p6|) with 
Tf(q,0) = B(q 2 ). When P^O then Tf(q,P) ^ B(q), and @ must be solved for Tf(q,P). 

Because the pion mass is small when m is small, we can perform an expansion of the P^ 
dependence in the kernel of (|7^). Since the analysis is Lorentz covariant we can, without loss of 
validity, choose to work in the rest frame with P = (0, im„), giving 

r(p) = \ml I J0^D(p-q)I(s)r(q)+ 



1 



A Swf D{p - q) 



rT(q) + ...., (71) 



' s(A(s) + m(s).e v (s)) 2 + (B{s) + m(s).e s ( S )) 2 
where 

100 = 6 (<J 2 V - 2(a s a' a + so v g' v ) - s(a s ^ - {a' s ) 2 ) - s 2 M ~ K,) 2 ) . (72) 



By using Fourier transforms the integral equation (71), now with explicit dependence on 



can be expressed in the form of a variational mass functional, 

MJT? = -- [**- ^ + 

T[ J U[T} 2 J (2tt)4 s(A(s) + m(s).e v (s)) 2 + (B(s) + m( S ).e s ( S )) 2 + 



6 r*r{*) 



2 



dx ^> (73) 



U[T] 2 J D(x) ' 
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which is equivalent to (j58| ) for the pion (but in the Feynman gauge), and in which 



/.[If = / ^I{s)T{qf. (74) 

The functional derivative SM^lT} 2 /ST(q) — reproduces (f7l|). The mass functional ( |73| ) and its 
minimisation is equivalent to the constituent pion BSE in the near chiral limit. To find an estimate 
for the minimum we need only note that the change in to 2 from its chiral limit value of zero will 
be of 1st order in to, while the change in T(q) from its chiral limit value B(q 2 ) will be of 2nd order 
in to. 

Hence to obtain to 2 to lowest order in to, we may replace T(q) by B(q 2 ) in (73), and we have 
that the constituent pion mass is given by 

a 48 f d 4 g e s (s)B(s) + se v (s)A(s) B(s) 2 



U[B} 2 J (2tt) 4 w sA(s) 2 +B(s) 2 sA(s) 2 + B( S y 



24 f S(S)2 ■ 6 f* x *§L + (Kn ?) (75) 



MB} 2 J (2tt) 4 sA(s) 2 + B(s) 2 U[B} 2 J D(x) 

However the pion mass has been shown to be zero in the chiral limit. This is confirmed as the two 
O(m ) terms in ( |7q ) cancel because of the identity ( p9|). Note that it might appear that f v would 
contribute an extra to dependence from its kernel in ([72|). However because the numerator in ([73]) 
is already of order to, this extra contribution must be of higher order in to. 

Hence we finally arrive at the analytic expression, to 0(m), for the constituent NG boson 
(mass) 2 from the solution of the BSE, 

2 48 f d 4 q e s (s)B(s) + se v (s)A(s) B{s) 2 



m * = MB?J Wr m[S) sA(s) 2 +B( S ) 2 sA( S ) 2 + B( S ) 2 + ° {m } ' (?6) 
It would appear that expression J7q ) is manifestly different to the conventional GMOR form in 

However here we generalise the identity found by Langfeld and Kettner 143] that shows these 
forms to be equivalent. Inserting ( |64| ) and ( |65| ) into (|6^), and expanding in powers of to(s), we 
obtain up to terms linear in m, and after using (p6[) with to = to eliminate the O(m ) terms, 



m(p 2 )e s (p 2 ) =m(p 2 ) + 4 —L-D(p - q) 



d 4 q „, , m(q 2 )e s (q 2 ) 



(2tt) 4 ^ HJ q 2 A(q 2 ) 2 + B(q 2 ) 2 
d 4 q n( , B(q 2 ) 2 2m(q 2 )e s (q 2 ) 



(2tt) 4 ^ ^' {q 2 A{q 2 ) 2 + B{q 2 ) 2 ) 2 

dS B(q 2 )A(q 2 )2m(q 2 )q 2 e v (q 2 ) 
-D(p-q)- 



(2tt) 4 ^ (q 2 A(q 2 ) 2 + B(q 2 ) 2 

(77) 

We now multiply ( |77j) t hroughout by B(p 2 )/(p 2 A(p 2 ) 2 + B(p 2 ) 2 ), and integrate wrt p. Using 
again the chiral limit of (|66|) there is some cancellation of terms, and we are left with a generalised 
Langfeld-Kettner identity 

dip B{p 2 ) 2 ( B(p 2 )m(p 2 )e s (p 2 ) , p 2 A(p 2 )m(p 2 )e v (p 2 ) 



p 2 A(p 2 ) 2 + B(p 2 ) 2 \p 2 A(p 2 ) 2 +B(p 2 ) 2 p 2 A{p 2 ) 2 + B{p 2 ) 2 

m(p 2 )B(p 2 
'p 2 A{p 2 ) 2 +B{p 2 ) 

Remarkably we see on using this identity that ((7^) is identical to ( |6l| ) 



d P -,2 A( n 2\2 , nf n 2\2- ( 78 ) 
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8 Constituent Nucleon 



We report here progress in calculating the constituent nucleon which emerges from the GCM hadro- 
nisation as a three constituent-quark state, bound by the effective gluon correlator. This nucleon 
state is treated as a correlation between a constituent quark and a constituent diquark subcorre- 
lation in the separable Faddeev approach. The first such computation was in 1989 |5(J and used 
a rank-1 description of the scalar diquark; this yielded a constituent core mass of approximately 
1.3GeV, which was very close to the expected core mass |5l[]. A full calculation of the nucleon core 
mass is particularly difficult and has yet to be attempted, as it requires the inclusion of the con- 
stituent quark correlators and the various constituent diquark correlators, but particularly that of 
the scalar + and vector 1 + diquarks. This nucleon core state then has its mass further reduced by 
pion dressing. As preparation for these extensive ab initio computations we have been monitoring 
the response of the nucleon core to the underlying low energy quark-gluon processes by computing 
the quark - scalar-diquark nucleon core state. The results are shown for GCM95 (rank 2), GCM97 
(rank 3) and GCM98 (rank 3) in Table 2. We now briefly outline the present procedures used in 



these studies. Working in a Euclidean metric the equation for the spin 2 
separable components (each a spinor) is 



i + 



nucleon form factor 



-E 



d q 



(27T) 



1 



2 -3a 



1 



PY)T 1 {{q+-p + 



3a 



Pf 



.G((2a -l)P-p- q)G((l - a)P + q)Z jk ((aP - q) 2 )^ k (q; P), 



(79) 



where the ^ k are defined in terms of an arbitrary momentum partitioning parameter a. Here the 
scalar diquark correlator is modelled using the separable form 



A(q,p,P) = '£r i (q)Z ij (P)T j (p), 



(80) 



where q and p are the relative quark momenta, and P is the diquark momentum. The matrix Zij 
is determined from the diquark version of the ladder BS correlator (^3|) when the gluon correlator 
separable expansion (^) is used, and only the 0(4) invariant part is retained. This diquark 
separable expansion is not to be confused with the formal expansion in (pl|). This nucleon core 
equation requires careful determination of its only ingredients, the quark correlator and the diquark 
correlator, and particularly its vertex functions. These are determined by solving the diquark BSE 
using the separable representation of the effective gluon correlator. No integration cutoffs are 
required. 

We seek solutions to ([79]) (which being a homogeneous linear equation will only have solutions 
for particular P 2 = —Mq) which give the nucleon core mass. We work in the rest frame of the 
nucleon and accordingly set P = (0, iM ). With the above choices of T, G, Z and P, ( f79|) enjoys a 
spatial 0(3) symmetry. A direct calculation shows that the integral operator commutes with the 
angular momentum operator J = L + S =i^xp + iir, so we take the ^ k to be one of the general 



L = 0, S 



i + 



states 



cf P 

IpI 



Uk(p) 
Vk(p) 



cr p 

IpI 



u k (p) 
Vkip) 



(81) 



where u k and v k are functions only of p4 and |p|. Equation (|79|) then becomes 



Viijp) 



-K ik (p A , |p| ; £? 4 , |q| ;p-q) 



Uk{q) 
Vk(q) 



(82) 
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Figure 11: The real and imaginary parts of the ?u(/>3,p4) and Vi(p3,p4) components (corresponding to 
the Fi in (BOh) of the nucleon wave function (scale: 10 units = 0.57GeV.) 
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where (G x ) 22 = [-((2o-l)M +t(p4+«4))A(«i) + Bfa)]^, 

G 2 : 



S2 A 2 (,S 2 )+B 2 ( S2 ) 

((l-a)M -%)A(s 2 ) + £(s 2 ) - |q| A(s 2 ) 

|q| A(s 2 ) -((l-a)M -%)A(s 2 ) + B(s 2 ) 

and the arguments of the quark and diquark correlators are 
si = {p± + qi - {2a - l)iM ) 2 + |p + q| 2 , 

s 2 = {qi + (1 - a)iM Q ) 2 + |q| 2 , } (83) 

s 3 = (g* - iaM ) 2 + |q| 2 . 



2G 



Equation ( p2|) reduces to n (the rank of the gluon correlator modelling) coupled two-dimensional 
integral equations after performing one trivial and one numerical angle integrations. We search 
for eigenvectors by introducing an eigenvalue A (Mo) and changing Mo until the eigenvalue A = 1. 
Values for Mo are in Table 2, and U\ (733,^4) and V\ (p3,Pi) are shown in Fig.[ll|. 

There have now been many quark-diquark Faddeev studies of the nucleon |32[ [53[ [54], [S5|, |5(| 
57, 58, E9I pol En], with applications to hadronic form factors |62, p3|. 



9 Pion Dressing of Constituent Nucleon 

The GCM hadronisation leads directly to the formalism for dressing the constituent nucleon by 



mesons as noted in sect. 3.5. The full determination of this nucleon state is not yet completed, 
however we illustrate here the nature of the meson dressing calculations. While chiral symmetry 
mandates the on-mass-shell pion- nucleon coupling, as in ( |39| ) , it is clearly necessary to include the 
ttNN form factor in calculating loop processes. The nucleon correlator Qn is then determined by 
the Euclidean metric DSE: 

g N \p)=^ + M + 3^^ J ^4 rp _ Q)2 + m 2 F ( P - Q> Q)i^GN(Q)il 5 To(Q - P, -Q), (84) 

where y 11 and j- are the core and dressed ttNN couplings, and Gn(P) has the form 

g N (P) = {iA N (P 2 )P.~, + B N (P 2 ) + M )-\ (85) 
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Figure 12: (a) Mass M of the dressed nucleon for various constituent nucleon masses Mo(GeV), for irNN 
form factor parameter /3 = 0.55GeV (longdash), /3 = 0.65GeV (solid) and /3 = 0.75GeV (shortdash), (b) 
Plots of the MIT bag constant B 1//4 (cr) for values of the scalar field < a < 1. The curves show results 
for GCM98-solid, GCM95-dashed and GCM97-shortdash. 

Here Mo is the constituent nucleon mass, while the mass of the dressed or physical nucleon M 
is given by the mass-shell condition P 2 An(P 2 ) 2 + £?tv(-P 2 ) 2 |p 2 =-m 2 = 0. This zero is situated in 
the time-like region and is determined by analytic continuation from the Euclidean support P 2 > 
in (Q). As M occurs in one of the couplings when solving ( |S4|) we must find self-consistency for 
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the value of M. This adds to the non-linearity of this DSE. We model the dressed and constituent 
form factors by the same separable approximation, 



r(p Q) = r Q (p, Q) = (86) 

(.1 + ^rJU + jz) 

in which P is the pion momentum and Q the nucleon momentum, where the parameters j3\ and 
/?2 are computable using the nucleon structure from sect.||, but here we show results, in Fig.|l2](a), 
for three typical values with f3\ = 02 ■ Despite its non-linearity ( [si] ) converges after a few itera- 
tions, indicating that the pion dressing of the nucleon core involves only a small number of pions. 
Nevertheless due to the low mass of the pion the nucleon mass shift of typically some 300MeV is 
significant, a result that is also seen in MIT bag modelling p4|. 



10 Connections to Other Models 



By further approximations the GCM provides a derivation of many other models, and so makes 
it possible to link these models to QCD as illustrated in Fig.|l|. A particular feature of this GCM 
linking is that it can predict the values of many of the phenomenological parameters occurring in 
these models, and relate their values to the underlying low-energy quark-gluon processes. Here we 
briefly indicate the connection to some of these models. 

The Nambu-Jona-Lasinio model (NJL) is an obvious special case of the GCM. Formally the 
NJL model jli} is the contact interaction limit of the GCM: D )iu {x — y) ~ + g 2 5 liv 5(x — y) or 
D(p) — > g 2 . But in the CQ equation the contact limit is undefined because it leads to divergences 
in ( |54| ) and (|55|). A cutoff A is then always introduced in NJL computations, which is equivalent 
to using the 'step-function' gluon correlator D(p) — > g 2 9(q 2 — A 2 ). Hence the NJL model is the 
GCM but with a box-shaped D(p), rather than a 'running' D(p). 

As we have already noted in sect. 3.4 the GCM provides a comprehensive derivation of the 
NG sector effective action. This is the ChPT effective action |55| , but with the added insight that 
all coefficients are given by explicit and convergent integrals in terms of A and B, which are in 
turn determined by D^. The higher order terms contribute to itit scattering, and the sensitivity 
of these to changes in P M „ are shown in Table 2. The GCM formalism also provides the non- 
local NG-baryon effective action, and this leads to finite values of observables, and so obviates the 
non-renormalisability problems that plague the local-ChPT phenomenology. 

While the GCM hadronisation is the main result, at an intermediate stage one obtains Jl| 
extended meson Quark-Meson Coupling type models (QMC) Applying mean field techniques 
to the GCM quark-meson coupling effective action leads to soliton type models, which have been 
studied in detail in j3lj and the significance of the extended mesons demonstrated. From the 
soliton models a further ansatz for the form of the soliton leads to the MIT and Cloudy Bag 
Model (CBM) ^37|. In |6£| baryons are modelled as hybrids of solitons and three quark bound 
states. An expression for the MIT Bag constant, derived from the GCM is jlj 



B = 



12tt 2 

(2*y 



sds 



ln( 



A 2 (s)s + B 2 {s) 
A^sja 



)- 



B 2 (s) 



A 2 (s)s + B 2 (s) 



(87) 



which is based on the energy density for complete restoration of the chirally symmetric pcrturbative 
configuration inside a cavity. This bag constant is for core states because no meson cloud effect is 
included. With a mean field description of the pion sector via <r(x), which describes the isoscalar 
part of a(x)V(x), where a{x) is a 'radial' field multiplying the NG boson field V(x) (see sect. 3.4)), 
B becomes, for constant a, 



B{a) 



12tt 2 
(2nY 



sds 



ln( 



A 2 {s)s + a 2 B 2 {s) 
A 2 (s)s 



a 2 B 2 {s) 



A 2 (s)s + B 2 (s) 
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which reduces to (87) when a = 1, being the non-perturbative field external to an isolated nucleon 
core. Here a < 1 describes a partial restoration of chiral symmetry outside of the core. Using 
the gluon correlator discussed in sect.^ we obtain the plot of 2? 1//4 (<t) shown in Fig.[l2](b), for the 
three GCM gluon correlators. Again we emphasize that the GCM provided not only the MIT 
bag phenomenology but also the value of the MIT bag constant. Dressing of the nucleon core by 
mesons, using a mean field modelling, is partly described by a reduction in a in the surface region, 
causing a reduction in the nucleon mass. 

However in nuclei a mean meson field description |]9| [7(J means that a is even further reduced 
outside of the nucleons, and the effective bag constant is further reduced. The a field can model 
in part correlated tttt exchanges and, along with the to meson field, is believed to be important 
in a mean field modelling of nuclei. In jnj it has been argued that the reduction of the effective 
bag constant for nucleons inside nuclei is essential to the recovery of features of relativistic nuclear 
phenomenology. The GCM thus allows B(a) and details of relativistic nuclear phenomenology to 
be directly related to the low energy quark-gluon processes that have been extracted from low 
energy meson data. 



11 Conclusion 

The GCM has turned out to be a very efficacious model of QCD when applied to low energy 
hadronic processes. This success appears to arise from a feature of QCD that might be thought to 
make low energy hadronic physics too difficult to sustain fundamental analytical models, namely 
the strength and number of gluonic processes in the IR regime. However their very strength seems 
to lead to an IR saturation effect in which the hadronic processes become somewhat insensitive to 
details of these gluonic processes. This fortuitous circumstance probably also explains why there 
are a considerable number of seemingly different but apparently successful hadronic models. The 
GCM appears to most successfully incorporate the manifestations of this simplifying feature of 
QCD. It does so by being itself a well-defined quantum field theory in which the consequences 
of the dynamical breaking of chiral symmetry are automatic and significant. It also supports the 
powerful hadronisation analysis from which the appropriate dynamical variables for low energy 
hadronic processes naturally emerge. A key part of this hadronisation is a bilocal meson-diquark 
bosonisation of the GCM. It is through this non-local bosonisation that we avoid the spurious 
introduction of a non-renormalisable effective action for the hadrons. Because of this we can 
uniquely relate the numerical values of numerous hadronic observables to the underlying modelling 
of the quark-gluon processes. This procedure is so robust that recent progress is already seeing 
the comparison of lattice-determined low-energy quark-gluon processes with those extracted from 
experimental data. Until now the GCM has mainly been applied to the meson sector, however 
as reported here work on an ab initio computation of the nucleon properties within the GCM is 
now well advanced. The nucleon is a complicated state to study not only because of its three 
quark character but also because the mesonic fluctuations play a significant role. The study of 
the baryonic sector of the GCM will provide a rich field of phenomena in which complex hadronic 
processes may be determined in the context of properly computable quantum amplitudes devoid 
of the non-renormalisability problems. 
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